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Abstract 



We explore some important implications of the KamLAND measurment on 
the lepton flavor mixing matrix V and the neutrino mass matrix M. The 
model-independent constraints on nine matrix elements of V are obtained 
to a reasonable degree of accuracy. We find that nine two-zero textures of 
M are compatible with current experimental data, but two of them are only 
marginally allowed. Instructive predictions are given for the absolute neutrino 
masses, Majorana phases of CP violation, effective masses of the tritium beta 
decay and neutrinoless double beta decay. 
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I. INTRODUCTION 



The KamLAND experiment |T| turns to confirm the large-mixing-angle (LMA) Mikheyev- 
Smirnov-Wolfenstein (MSW) solution to the long-standing solar neutrino problem. In 
addition, the K2K long-baseline experiment has unambiguously observed a reduction of 

flux and a distortion of the energy spectrum. These new measurements, together with the 
compelling SNO evidence [^] for the flavor conversion of solar z/g neutrinos and the Super- 
Kamiokande evidecne ||^ for the deficit of atmospheric neutrinos, convinces us that the 
hypothesis of neutrino oscillations is actually correct! We are then led to the conclusion that 
neutrinos are massive and lepton flavors are mixed. 

The mixing of lepton flavors means a mismatch between neutrino mass eigenstates 
(i^i, i/2! i^s) and neutrino flavor eigenstates {i'f.,i'^,i'r) in the basis where the charged lep- 
ton mass matrix is diagonal: 
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Taking account of the unitarity of V, one may reversely express |Vei|, |V^2|, iKsl, 
iKsl in terms of 6'sun, O^tm and 6'chz 0: 
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Current experimental information on 6'sun, ^atm and ^chz allows us to get very instructive 
constraints on the lepton flavor mixing matrix V. One purpose of this paper is therefore to 
examine how accurately we can recast V from the present KamLAND, K2K, SNO, Super- 
Kamiokande and CHOOZ measurements. 

Another purpose of this paper is to confront two-zero textures of the neutrino mass 
matrix with the new KamLAND data, so as to single out the most favorable texture(s) in 
phenomenology. In the flavor basis chosen above, the Majorana neutrino mass matrix can 
be written as 
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where m^ (for i = 1,2,3) are physical masses of three neutrinos. In the assumption of the 
LMA solution for solar neutrino oscillations, a classification of M with two vanishing entries 
has been done in an analytically approximate way P,^. Our present analysis is different 
from the previous ones in two important aspects: (1) we carry out a careful numerical 
analysis of every two-zero pattern of the neutrino mass matrix M to pin down its complete 
parameter space, because simple analytical approximations are sometimes unable to reveal 
the whole regions of relevant parameters allowed by new experimental data; (2) we present 
the quantitative predictions for allowed ranges of the absolute neutrino masses, the Majorana 
phases of CP violation, and the effective masses of the tritium beta decay ((m)e) and 
neutrinoless double beta decay ((m)ee). 

The remaining part of this paper is organized as follows. With the help of current data, 
we derive the model-independent constraints on nine elements of the lepton flavor mixing 
matrix in section II. Section III is devoted to a detailed analysis of the parameter space for 
every two-zero texture of the neutrino mass matrix. We obtain instructive predictions for 
the neutrino mass spectrum and Majorana phases of CP violation as well as (m)e and (m)ee 
in section IV. Finally a brief summary is given in section V. 



II. CONSTRAINTS ON THE LEPTON FLAVOR MIXING MATRIX 



As already shown in Eq. (3), five matrix elements of V can be determined or constrained 
from current experimental data. The other four matrix elements (|V^i|, |V^2|; and 
\Vr2\) are entirely unrestricted, however, unless one of them or the rephasing invariant of 



CP violation of V (defined as J [|T0[) is measured. A realistic way to get rough but useful 
constraints on those four unknown elements is to allow the Dirac phase of CP violation 
in V to vary between and vr |jll|], such that one can find out the maximal and minimal 
magnitudes of each matrix element. To see this point more clearly, we adopt the standard 
parametrization V = UP [O, where 
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with Sr = sin 6r, Cr = cos^r, and so on; and 
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The advantage of this representation is that the neutrinoless double beta decay is associated 
with the Majorana phases p and a, while CP violation in normal neutrino oscillations 
depends separately on the Dirac phase 5. Note that three mixing angles {0x,9y,9z), which 
are all arranged to lie in the first quadrant, can be written as 
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It is then straightforward to obtain 
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Varying the Dirac phase 5 from to tt, we are led to the most generous ranges of |V^i|, | V^2|, 
iKil and IK2I: 
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Note that the lower and upper bounds of each matrix element turn to coincide with each 
other in the limit IK3I — > 0. Because of the smallness of iKsl, the ranges obtained in Eq. (9) 
should be quite restrictive. Hence it makes sense to recast the lepton flavor mixing matrix 
even in the absence of any experimental information on CP violation. 

In view of the present experimental data from KamLAND [|l| , K2K ^ , SNO , Super- 
Kamiokande § and CHOOZ H, we have 0.25 < sin^ O^nn < 0.40 H, 0.92 < sin^ 2^atm < 1-0, 
and < sin^ 26'chz < 0.1 at the 90% confidence level. Namely, 



30.0° < ^sun < 39.2° , 
36.8° < 0atm < 53.2° 
0° < ^^chz < 9.2° . 



(10) 



Using these inputs, we calculate the numerical ranges of |Ki|, IK2I, IK3I) I^AisI and IK-3I 
from Eq. (3). Then the allowed ranges of |V^i|, |V^2|, and IK2I can be found with the 
help of Eq. (9). Our numerical results are summarized as 
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This result is certainly more restrictive than that obtained in Ref. ||TT| before the KamLAND 
measurement. 

Note that the rephasing invariant of CP violation reads as follows: 
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The term proportional to 4 sin^ 6'atm sin^ 6'chz in J, which may correct the leading term up 
to 5% (for 6'atm = 45° and 6'chz = 9°), were not taken into account in Ref. |Tl|]. By use 



of Eq. (10), we find J < 0.039 sin 5. This result implies that the magnitude of J can 
maximally be 0.039, leading probably to observable CP-violating effects in long-baseline 
neutrino oscillations. 



III. TWO-ZERO TEXTURES OF THE NEUTRINO MASS MATRIX 

The symmetric neutrino mass matrix M totally has six independent complex entries. If 
two of them vanish, i.e., Mab = Mpq = 0, we obtain two constraint equations: 

mif/pit/.ie^^^ + m2Up2U,2e^"' + m^U,,^U,^ = , (13) 
where a, 6, p and q run over e, /i and r, but (p, q) ^ (a, h). Solving Eq. (12), we arrive at 



'^^2ip _ UazUbzUp2Uq2 - Ua2Ub2Up?,UqZ 

1^3 Ua2Ub2UpiUgi - UalUblUp2Uq2 ' 

"^2 ^2icT _ UalUblUpzUq3 — UazUblJJplUql /-j^^n 

^3 Ua2Ub2UpiUqi - UalUblUp2Uq2 ' 

This result implies that two neutrino mass ratios (mi/ms, 777,2/^3) and two Majorana-type 
CP-violating phases (p, a) can fully be determined in terms of three mixing angles {dx, 6y, 0^) 
and the Dirac-type CP- violating phase (S). Thus one may examine whether a two- zero 
texture of M is empirically acceptable or not by comparing its prediction for the ratio of 
two neutrino mass-squared differences with the result extracted from current experimental 
data on solar and atmospheric neutrino oscillations: 



P. - ^-Pf^ - ^ . (15) 
\mi-mi\ Ami,^ 



Considering the LMA MSW solution confirmed by the KamLAND measurement, we have 
5.9 X 10"^ eV^ < Am^^^ < 8.8 x 10"^ eV^ [|l|,|ig at the 90% confidence level. In addition, 



5 



we have 1.6 x 10 ^ eV^ < Aml^^ < 3.9 x 10 ^ eV^ at the 90% confidence level Thus we 



arrive at 1.5 x 10~^ < -Rj^ < 5.5 x 10~^. The allowed ranges of three mixing angles 6^ ^ dsnn, 
~ ^atm and 9z ~ 6'chz have been given in Eq. (10). There is no experimental constraint 
on the CP-violating phase 5. Hence we simply take 5 from 0° to 360° in our numerical 
calculations. 

There are totally fifteen distinct topologies for the structure of M with two independent 
vanishing entries, as shown in Tables 1 and 2. We work out the explicit expressions of 
{mi/m^)e^'^'' and (m2/m3)e^*°" for each pattern of M by use of Eq. (14), and list the results 
in the same tables |]TB|]. With the input values of O^, Oy, 9^ and 6 mentioned above, we 
calculate the ratio R^, and examine whether it is in the range allowed by current data. This 
criterion has been used in Refs. to pick the phenomenologically favored patterns of M 
in the LMA case. 

Nine of the fifteen two-zero textures of M listed in Table 1 are found to be in accord 
with the LMA solution as well as the atmospheric neutrino data. They can be classified 
into four categories: A (with Ai and A2), B (with Bi, B2, B3 and B4), C and D (with 
Di and D2). The point of this classification is that the textures of M in each category 
result in similar physical consequences, which are almost indistinguishable in practice. The 
other six patterns of M (categories E and F) listed in Table 2 cannot coincide with current 
experimental data. In particular, the exact neutrino mass degeneracy (mi = m2 = ms) is 
predicted from three textures of M belonging to category F. 

Now let us focus on patterns Ai, Bi, C and Di as four typical examples for numerical 
illustration. Our results for sin^ 20chz versus 6 and 9y versus 9x are shown Figs. 1-4. Some 
comments are in order. 

(1) For pattern Ai, arbitrary values of 6 are allowed if sin^ 26'chz is large enough (> 0.014). 
The mixing angles 9^ and 9y may take any values in the ranges allowed by current data. 
Therefore we conclude that pattern Ai is favored in phenomenology with little fine-tuning. 
A similar conclusion can be drawn for pattern A2. 

(2) For pattern Bi, 5 is essentially unconstrained if sin^ 26'chz is extremely close to zero; 
and only 6 around 90° or 270° is acceptable if sin^ 26'chz deviates somehow from zero. Except 
9y 7^ 45°, there is no further constraint on the parameter space of {9^, 9y). We conclude that 
pattern Bi with maximal CP violation (i.e., sin 5 ~ ±1) is phenomenologically favored. So 
are patterns B2, B3 and B4. 

(3) For pattern C, 5 = 90° or 5 = 270° is forbidden. Furthermore, 9y = 45° is forbidden. 
We see that the allowed parameter space of {S, 9chz) and that of {9^, 9y) are rather large. 
Hence pattern C is also favored in phenomenology. 

(4) For pattern Di, S is restricted to be around 0° or 360°. In particular, the region 
90° < 5 < 270° is entirely excluded, sin^ 2^chz > 0.084 holds for the allowed range of 
6. Different from patterns Ai, Bi and C, pattern Di requires relatively strong correlation 
between 9^ and 9y (e.g., small values of 9y are associated with large values of 9^ in the allowed 
parameter space). In this sense, we argue that pattern Di is less natural in phenomenology, 
although it has not been ruled out by current experimental data. A similar argument can 
be made for pattern D2. 

It is worth remarking that patterns Di and D2 were not included into the phenomeno- 
logically allowed patterns of M in the previous classification PJ^, where only analytical 
approximations were made. Our numerical analysis shows that these two patterns are 
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marginally allowed by current data. Here we have also explored some interesting details 
of the parameter space for every favorable texture, which could not be seen from simple 
analytical approximations. 



IV. NUMERICAL PREDICTIONS AND FURTHER DISCUSSIONS 

A two- zero texture of M has a number of interesting predictions, in particular, for the 
absolute neutrino masses and the Majorana phases of CP violation 0. With the help of 
Eq. (14), one may calculate the mass ratios mi/m^ and 7712/7713 as well as the Majorana 
phases p and a. The absolute neutrino mass can be determined from 
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(16) 



Therefore a full determination of the mass spectrum of three neutrinos is actually possible. 
Then we may obtain definite predictions for the effective mass of the tritium beta decay. 



(777), 



2 2 I 2 2 I 2 
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and that of the neutrinoless double beta decay. 
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It is clear that the Dirac phase 5 has no contribution to {m)ee- Note that CP- and T- 
violating asymmetries in normal neutrino oscillations are controlled by 6 or the rephasing- 



invariant parameter J 



sin (5. Whether (777)e and (777)^ 



can be measured 



remains an open question. The present experimental upper bounds are (77i)e < 2.2 eV 
and (777)ee < 0.35 eV |18] at the 90% confidence level. The proposed KATRIN experiment 
is possible to reach the sensitivity (77i)e ~ 0.3 eV []TU[, and a number of next-generation 
experiments for the neutrinoless double beta decay pOf is possible to probe (m)ee at the 
level of 10 meV to 50 meV. 

We perform a numerical calculation of m^jvn'^ versus mi/m^, o versus p, (?77)ee versus 
(77i)e, and J versus 7773 for patterns Ai, Bi, C and Di. The results are shown in Figs. 1-4. 
Some discussions are in order. 

(1) For pattern Ai, p ^ 6/2 or p ^ 6/2 — 180° and cr ^ p ± 90° hold in most cases. Two 
neutrino mass ratios lie in the ranges 0.033 < nii/m^ < 0.19 and 0.13 < 7772/7713 < 0.28, 
and the absolute value of 7773 is in the range 0.04 eV < 7773 < 0.065 eV. As (777)ee = 
is a direct consequence of texture Ai, we calculate the sum of three neutrino masses 
instead of {m)ee- The result is 0.047 eV < < 0.093 eV, in contrast with 0.003 eV < 
(777)e < 0.014 eV. The rephasing invariant of CP violation J is found to lie in the range 
—0.037 < J < 0.038. Similar predictions are expected for pattern A2. 

(2) For pattern Bi, p ^ a ^ 6 — 90° or p ^ a ~ 5 — 270° holds in most cases. Two 
neutrino mass ratios mi/m^ may lie either in the range 0.53 < rrti/m^ < 0.99 or in the 
range 1.01 < rrii/m^ < 1.88, and 7772/7773 may lie either in the range 0.53 < 7772/7713 < 0.99 
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or in the range 1.01 < 1712/1713 < 1.88. The value of ms is found to be in the range 
0.026 eV < m3 < 0.35 eV. Furthermore, we arrive at 0.027 eV < (m)e ~ (?7i)ee < 0.35 eV 
as well as —0.037 < J < 0.037. Similar results can be obtained for patterns B2, B3 and B4. 

(3) For pattern C, a ~ p when 6y approaches 45°; and there is no clear correlation 
between p and a for other values of 6y. Two neutrino mass ratios nii/m^ may lie either in 
the range 0.95 < mi/m^ < 0.99 or in the range 1.01 < nii/m^ < 5.4, and m2/m3 may lie 
either in the range 0.95 < m2/m3 < 0.99 or in the range 1.01 < 1712/1173 < 5.3. The value of 

is found to lie in the range 0.009 eV < < 0.35 eV. It is remarkable that {m)ee ~ '"1.3 
holds to a good degree of accuracy in the allowed space of those input parameters. We also 
obtain 0.04 eV < (m)^ < 0.35 eV and -0.037 < J < 0.037. 

(4) For pattern Di, p ^ 6 — 90° or p ^ 5 — 270° and cr ^ p ± 90° hold. Two neutrino 
mass ratios lie in the ranges 7.5 < nii/m^ < 8.8 and 7.35 < m2/m3 < 8.6, and the 
absolute value of is in the range 0.005 eV < < 0.008 eV. As for the tritium beta 
decay and neutrinoless double beta decay, we obtain 0.04 eV < (m)e < 0.062 eV and 
0.008 eV < (m)ee < 0.014 eV. The range of J is found to be -0.014 < J < 0.011. Similar 
predictions can straightforwardly be made for pattern D2. 

We see that there is no hope to measure both (m)e and {in)ee, if the neutrino mass 
matrix M takes pattern Ai or A2. As for categories B and C of M, the upper limit of (m)e 



is close to the sensitivity of the KATRIN experiment (~ 0.3 eV [1£]), and that of (m)ee is 
just below the current experimental bound 



V. SUMMARY 

In summary, we have discussed some implications of the KamLAND measurment on the 
lepton flavor mixing matrix V and the neutrino mass matrix M. The model-independent 
constraints on nine elements of V have been obtained up to a reasonable degree of accuracy. 
Nine two-zero textures of M are found to be compatible with current experimental data, 
but two of them are only marginally allowed. Instructive consequences of these phenomeno- 
logically favored textures of M on the absolute neutrino masses, Majorana phases of CP 
violation, (m)e and (m)ee are numerically explored. Our results will be very useful for model 



building |21|, in order to understand why neutrino masses are so tiny and why two of the 
lepton flavor mixing angles are so large. 

Finally it is worth remarking that a specific texture of lepton mass matrices may not be 
preserved to all orders or at any energy scales in the unspecified interactions from which 
lepton masses are generated. Nevertheless, those phenomenologically favored textures at low 
energy scales, no matter whether they are of the two-zero form or other forms, are possible 
to provide enlightening hints at the underlying dynamics of lepton mass generation at high 
energy scales. 

We would like to thank Y.F. Wang for stimulating discussions about the KamLAND 
measurement. One of us (Z.Z.X.) is also grateful to IPPP in University of Durham, where 
the paper was finalized, for its warm hospitality. This work was supported in part by 
National Natural Science Foundation of China. 
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TABLES 

TABLE I. Nine patterns of the neutrino mass matrix M with two independent vanishing entries, 
which arc com,patihle with the LMA solution and other empirical hypotheses. The analytical results 
for two ratios of three neutrino mass eigenvalues {m\/m^)e^^^ and {1712/ ms)e'^^'^ are given in terms 
of four flavor mixing parameters 9j:, 9y, 9^ and S. 
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m2g2i(T _ I SxCz _ Sx[sy—c'^)—2cxSyCySze^^ 



"^3 2sxCxSyCy — {s'^—c%)(^s^—c'^')sze''^+2sxCxSyCys'^e- 



2x5 



rrn^2ip _ _£z . C3:Sy 

ma Sa^Cy+CajSySze'* 

m2g2iiT _ . gi<5 

ma Sz CxCy—SxSySze*" 



mi g2jp _ _|_£| . 

77^3 ^2 SxSy CxCySz^^ 
2 

mS Sz CxSy+SxCySzC^^ 
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TABLE II. Six patterns of the neutrino mass matrix M with two independent vanishing entries, 
which arc incompatible with the LMA sohition and other empirical hypotheses. The analytical 
results for two ratios of three neutrino mass eigenvalues (mi/msje^*'' and (777-2/^3)6^*'^ are given 
in terms of four flavor mixing parameters 9x, 9y, 9z and S. 

Pattern of M Results of {mi/mz)e^'^P and (7712/7773)6^*'^ 



El : 



Eo 



Eq 






X 


X 


X 





X 


X 


X 


X 





X 


X 


X 


X 


X 


X 


X 








X 


X 


X 


X 





X 





X 


X 











X 


X 





X 


X 


X 





X 





X 





X 





X 


X 


X 





X 


X 











X 



m\ „2io 1 ^x^y{'^^ *z) CxCyS^(^CxCy 2sxSySze'-^^e oiS 

— 6 ^--^ « 

rmjiia _ J 1_ . c^s^[^c'^-s'z)-SxCyS'z\^SxCy+-zca:SySze- )e-'- ^2iS 

ma 



'1 • I " \ mL2ia 





„2 2 ( 2_ 


(si- 


^'x^^y~^^^^^^ ^y ^ z^^^ 

-SxCyS^ ^SxCy-\-'2CxSySz C^^^G "^^"^ 






G^^Cy-\-2,SxCxSySzG^ 


1 




^^)- 


-CxSyS ^(^CxSy-\-'2iSxCyS zG 


+ 7^ 


^x^y V z 


(si- 

si)- 

(c2 


~C^)Sy '2s xCxCy S zG^^ 

- SxSyS (^SxSy 2CxCyS zG ^^G ^ 



m\ Jlip 1 s^Cy\Cz—Sz)—CxSySz\CxSy+ASxCySze je oiS 

—6 y--^ e 

1712 2ia I 1 ^x'^y{^z~Sz )~SxSySz [^SxSy—zcxCySze'" Je 2iS 



mg^iP = - 1 _ slsyCy{ci-sl)+Cxsl[cxSyCy-Sx{sl-cl)sze'^]e 2i<5 ^2j5 
ms {cl-sl)syCy+SxCx{cl-sl)szei^ 

m2fp.ia ^ I 1 . <^lsycy{cl-sl)+sxsl[sxsycy+cx{sl-cl)szei^]e-'^i^ ^^iS 
ms (cl-sl)syCy+SxCx{cl-sl)sze'^ 



mi g2jp _ 
^2ia ^ I 



rni^2ip _ SxCy+CxSySze ^2i5 
ms SxCy+CxSySze+'''' 



Tn2 2ia CxCy—SxSySze ''^ „2iS 

SxSyS zG 



mi g2«P = SxSy-CxCySze '° ^21.5 
ms SxSy—CxCySze+'''' 

m2g2iiT _ CxSy+SxCySzC g2*i5 
ms CxSy+SxCySze+'^'' 
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